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Abstract. We present results based on the analysis of optical 
spectra of KUG 1031+398, a Narrow Line Seyfert 1 (NLS1) 
galaxy for which Mason et al. (1996) reported evidence for a 
line-emitting region "intermediate" (both in terms of velocity 
and density) between the conventional broad and narrow line 
regions (BLR and NLR, respectively). From our observations 
and modelling of the spectra, we get a consistent decomposi- 
tion of the line profiles into four components: an extended H II 
region with unresolved lines, two distinct Seyfert-type clouds 
identified with the NLR, and a relatively narrow "broad line" 
component emitting only Balmer lines but no forbidden lines. 
Therefore, and although we find this object to be exceptional in 
having line-emission from the BLR with almost the same width 
as the narrow lines, our interpretation of the data does not sup- 
port the existence of an "intermediate" line region (ILR). 

Key words: galaxies: active - galaxies: Seyfert - galaxies: in- 
dividual: KUG 1031+398 



1. Introduction 

1.1. The Intermediate Line Region 

It is commonly accepted that line-emission in AGNs comes 
from two well separated regions: one, compact (smaller than 
1 pc) and lying close to the central engine, has a high electron 
density (N c > 10 8 cm -3 ) and is responsible for the production 
of broad (FWHM ~ thousands of kms -1 ) permitted lines - 
the BLR; the other, more extended and lying further away from 
the central source (10-1000 pc), has lower electron densities 
(10 3 < iV < 10 5 cm -3 ) and emits lines with a lower velocity 
dispersion (~ hundreds of km s _1 ) - the NLR. 

A line-emission "gap" is usually observed between the two 
regions, most objects showing an optical spectrum which can 
be fitted by line profiles corresponding to clouds belonging to 
one or the other line-emitting regions. This line-emission gap 
can be explained by the presence of dust mixed with the gas 
(Netzer & Laor 1993). Nevertheless, the existence of an inter- 
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mediate region, both in terms of velocity and density, is ex- 
pected; in such a region, the [O III] lines would be partially col- 
lisionaly de-excitedQ and show substantially broadened wings 
(Shields 1978). This ILR should not be confused with the ILR 
found in QSOs by Brotherton et al. (1994), which is much 
smaller and denser, with a velocity dispersion of the order of 
2000 km s" 1 and density - 10 10 cm" 3 . 

Mason et al. (1996) presented high-resolution (2 A) optical 
spectroscopic observations of KUG 103 1+398. The model they 
used to fit the data revealed a line-emitting region with lines 
of intermediate width (FWHM — 1 000 kms -1 ); according to 
Mason et al., this region would dominate the Balmer lines pro- 
file, being also a significant contributor to the [O lll]AA4959, 
5007 lines, with a flux ratio A5007/H/3 = 1.4, suggesting an 
intermediate density. 

Osterbrock (1978) thought that he had detected, in a few 
Seyfert 1 galaxies, faint wings to the [O III] lines with essen- 
tially the same widths as the Balmer lines. Crenshaw & Peter- 
son (1986) and van Groningen & de Bruyn (1989) have found 
broad wings in the [O III] lines of a number of Seyfert 1 galax- 
ies, implying the presence, in these objects, of an ILR with a 
density of a few times 10 6 cm~ 3 , similar to the one reported in 
KUG 1031+398; however, all these objects show strong Fell 
emission, and the observed broad [O III] components could be 
due to an inaccurate removal of the Fe II blends (Boroson & 
Green 1992). 

In summary, although the presence in Seyfert 1 galaxies 
of emitting clouds with density intermediate between those of 
the "broad" and "narrow" components is not unexpected, no 
uncontroversial report of the existence of such intermediate 
components has ever been made to the best of our knowledge. 
Therefore, the claims by Mason et al. (1996) that the NLS1 
KUG 1031+398 shows evidence for an ILR induced us to con- 
duct new spectroscopic observations and modelling of its emis- 
sion-line features. 



1 If we make the assumptions that the excitation conditions in both 
the NLR and the ILR are the same and that in the ILR, A5007/H/3 
~ 1, the [O III]A5007 line is collisionaly de-excited by about a factor 
10. According to the formula given by Seaton (1975), this implies a 
density of 1-3 10 6 cm~ 3 if the electronic temperature is in the range 
1-3 10 4 K. 
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1.2. KUG 1031+398 

2RE J 1034+396 was found in the ROSAT Wide Field Camera 
all-sky extreme-ultraviolet survey (Pounds et al. 1993; Pye et 
al. 1995). It was identified by Shara et al. (1993) and Mason et 
al. (1995) with the compact UV-excess 15.0 mag galaxy KUG 
1031+398 (Takase & Miyauchi-Isobe 1987) at z = 0.042. This 
object has an intense soft X-ray emission with an unusually 
steep 2-10 keV power law of photon index T ~ 2.6 ±0.1 
(Pounds et al. 1995) and an even steeper 0.1-2.4 keV power 
law with r = 3.4 ± 0.3 (Puchnarewicz et al. 1995) or V = 4.4 
±0.1 (Rodriguez-Pascual et al. 1997). 

UV spectroscopy with the Hubble Space Telescope shows 
the Ly a profile to be complex, with a narrow (400 kms -1 
FWHM) and a broad (1600 kms" 1 FWHM) component 
(Puchnarewicz et al. 1998). The optical continuum is not po- 
larized (Breeveld & Puchnarewicz 1998). 

The broad component of the Balmer lines is relatively nar- 
row (FWHM ~ 1500 kms -1 ) and, consequently, this object 
has been classified as a NLS1 by Puchnarewicz et al. (1995). 
Narrow line Seyfert 1 galaxies are defined as Seyfert Is having 
"broad" Balmer lines narrower than 2 000 km s _1 FWHM (Os- 
terbrock 1987). Most NLSls have a steep soft (< 1 keV) X-ray 
component and, conversely, most ultra-soft X-ray sources are 
associated with a NLS1 (Puchnarewicz et al. 1992; Greiner et 
al. 1996; Boiler et al. 1996; Wang et al. 1996). 

2. Observations 

Spectroscopic observations of KUG 1031+398 were carried 
out with the spectrograph CARELEC (Lemaitre et al. 1989) 
attached to the Cassegrain focus of the Observatoire de 
Haute Provence (OHP) 1.93 m telescope. The detector was a 
512x512 pixels, 27x27 pm Tektronix CCD. We used a 600 
1 mm -1 grating giving a dispersion of 66 A mm -1 . On January 
10, 1997 we obtained a 20 min exposure in the range A A 6175- 
7075 A, on March 4, a 20 min exposure in the range AA 4780— 
5780 A, and three more on March 5. 

The slit width was 2'.'1, corresponding to a projected slit 
width on the detector of 52 pm, or 1.9 pixel. The resolution, as 
measured on the night sky emission lines, was 3.4 A FWHM in 
the blue and 3.5 A in the red. In both cases the galaxy nucleus 
was centered on the slit and 3 columns of the CCD (~ 3'.'2) 
were extracted, corresponding to ~ 4 kpc at the distance of the 
galaxy (with H Q = 50 kms -1 Mpc -1 ). 

The spectra were flux calibrated using the standard stars EG 
247 (Oke 1974) and Feige 66 (Massey et al. 1988), observed 
with the same instrumental settings; these standards were also 
used to correct the red spectrum for the atmospheric B band at 
A6867 A (Fig. |]b). 

3. Analysis 

3.1. Methodology 

Positive correlations between line-widths and ionization poten- 
tials/critical densities have been observed in the narrow line re- 



gion of many Seyfert galaxies. Negative correlations are also 
sometimes found. A positive correlation implies that the den- 
sity and/or ionization parameter gradually increases inward in 
the NLR of these objects. In the case of correlations with criti- 
cal densities, the observed values range from ~ 10 3 to ~ 3 10 6 
cm~ 3 (Filippenko & Sargent 1985). 

Sub-structures were found in the narrow line profiles of 
most objects suggesting that the line emitting region is a col- 
lection of individual clouds in motion relative to each other and 
producing different parts of the line profiles (see for instance 
Veilleux 1991, Espey et al. 1994 or Ferguson et al. 1997). 

These findings induced us to assume that each of these 
clouds is characterized by a single density and that lines com- 
ing from the same emission-region should have the same pro- 
file and mean velocity. We therefore tried to model the spectra 
of KUG 1031+398 with the smallest possible number of line 
sets, each set including three Gaussians (modelling Ha and the 
[Nn]AA6548, 6583 lines, or H/3 and the [Om]AA4959, 5007 
lines) having the same velocity shift and width, with the ad- 
ditional constraint that the intensity ratio of the two [N II] (re- 
spectively [O III]) lines was taken to be equal to the theoret- 
ical value of 3.00 (respectively 2.96) (Osterbrock 1974). In a 
physically meaningful and self-consistent model, the compo- 
nents found when fitting the blue and red spectra should have 
velocity shifts and widths compatible within the measurement 
errors. 

3.2. The narrow line-core components 

The spectra were deredshifted assuming z = 0.0434 (Fig. 0a 
and b) and analysed in terms of Gaussian components as de- 
scribed above. We discovered first that the core of the lines 
could not be fitted by a single set of narrow Gaussian pro- 
files. To get a satisfactory fit, two sets of Gaussian compo- 
nents are needed: the first, unresolved (and subsequently taken 
as the origin of the velocity scales) has A6583/Ha = 0.55 and 
A5007 /H/3 = 1 .27 and corresponds to a H II region (Fig. 0); the 
second is resolved (FWHM ~ 350 kms -1 , corrected for the 
instrumental broadening), blueshifted by ~ 95 km s _1 with re- 
spect to the narrower components and has line intensity ratios 
typical of a Seyfert 2 (A6583 /Ha = 0.84, A5007 /H/3 = 10.2). 

3.3. The "broad" components 

At this stage, we removed from the blue and red spectra the 
best fitting line-core (the H II region and the Seyfert 2 nebulos- 
ity, Fig. 0: and d), obtaining two spectra we shall call "original 
data minus core". The blue one was then fitted with a broad H/3 
Gaussian component and two sets of three components mod- 
elling the narrow H/3 and [O m]AA4959, 5007 lines. The re- 
sult is very suggestive: one set has a strong H/3 line and very 
weak negative [O III] components, while the other set displays 
a strong [O III] contribution and a weak negative H/3 compo- 
nent, showing that we have in fact a H/3 component with no 
associated [O III] emission and [O III] lines with a very weak 
(undetected) associated H/3; in other words, the region produc- 
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Fig.l. Blue (a) and red (b) 

spectra of KUG 1031+398 in 
the rest frame; in (b) we also 
give the spectrum before cor- 
recting for the atmospheric ab- 
sorption (dotted line). The nar- 
row line-core components (c 
and d) were fitted with Gaus- 
sians and subtracted from the 
original data, the result being 
shown in (e) to (h). In (e) 
and (f), we show our best fit 
(solid line) together with the 
data points (crosses); the lower 
solid lines represent the residu- 
als. 



ing the H/3 line does not emit forbidden lines, while the [O III] 
emitting region has a high A5007/H/3 ratio, which are, respec- 
tively, the characteristics of the "broad" and "narrow" line re- 
gions in Seyfert 1 galaxies. 

Having these results in mind, we optimized this last fit by 
using a Lorentzian profile for the H/3 line, with no associated 
[O III] emission, and a set of three Gaussians for the remaining 
contribution coming from the "narrow" components (this is not 
the first time Lorentzian profiles are used to fit AGN emission 
lines; for example Stiiwe et al. (1992) found that, in the case of 
NGC 4258, the narrow lines were better fitted by Lorentzians, 
rather than Gaussians). 

The best fit is presented in Fig. |l^: in this model, the 
flux of the "narrow" (Gaussian) H/3 component is only 9% of 
the "broad" (Lorentzian) H/3 component. The H/3 Lorentzian 
component is blueshifted by 150 kms -1 with a width of 915 
kms -1 . The Gaussian components are blueshifted by ~ 395 
km s ~ 1 and their width is ~ 1 1 1 5 km s ~ 1 1 . Mason et al. found 
a FWHM = 1 030 ±150 kms -1 for this component which is 



blueshifted by 240 ± 30 kms -1 . This blueshift, however, is 
measured with respect to the [O III] lines core which is dom- 
inated by the Seyfert 2 cloud, itself blueshifted by 95 kms -1 
with respect to the H II region; the blueshift of Mason et al.'s 
intermediate component is, therefore, 240 + 95 = 335 kms -1 , 
in agreement with our value of 395 km s ~ 1 . 

In Seyfert 1 galaxies, the Balmer decrement of the broad 
component is never smaller than that of the narrow component; 
in the present case, we therefore expect the "narrow" Ha com- 
ponent flux to be less than 9% of the "broad" Ha component 
flux. Moreover, the [Nll]A6583 line flux is, in Seyfert galax- 
ies, equal or smaller than the narrow Ha component flux. So, 
in KUG 103 1+398, we expect the narrow lines to be quite weak 
compared to the "broad" Ha component, and we fitted the 
"original data minus core" red spectrum with a single Loren- 
tzian profile of 1 205 kms -1 FWHM, blueshifted by 65 km s -1 
with respect to the H II region. This fit is shown in Fig. |l]f. 

Another model allowing, in addition, for a set of Ha and 
[N II] Gaussian components was also tested, resulting in a fit 
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of similar quality; the very small A6583 /Ha ratio observed for 
this solution (0.2), shows that the nitrogen lines may be consid- 
ered as undetectable. 

4. Results and Discussion 

Our new observations and modelling of KUG 103 1+398 yield a 
consistent decomposition of the emission-line profile into four 
components (see Table [[]): 

1. An extended H II region with unresolved lines (Fig. ||); 

2. A first Seyfert-type cloud with relatively narrow lines (~ 
350 kms -1 FWHM), blueshifted by 95 kmr 1 , belonging 
to the NLR; 

3. A second Seyfert-type component with somewhat broader 
lines, blueshifted by ~ 395 kms -1 ; the width of the lines 
in this component (~ 1 115 kms -1 FWHM), which may 
seem large for a Seyfert 2, is not exceptional as the FWHM 
of the lines in the prototype Seyfert 2 galaxy NGC 1068 is 
- 1 670 kms" 1 (Alloin et al. 1983). Only the [Olll] lines 
are observed in this component, with A5007/H/3 ~ 6.1, a 
line ratio characteristic of NLRs; 

4. Finally, a Narrow Line Seyfert 1 cloud with lines well fit- 
ted by a Lorentzian profile of ~ 1060 kms -1 FWHM, 
blueshifted by 105 kms -1 . 

Our analysis shows that the emission line spectrum of KUG 
1031+398 can be satisfactorily decomposed in a set of compo- 
nents which have line ratios characteristics of H II regions or 
conventional NLR or BLR clouds, whithout the need to invoke 
the presence of an ILR characterized by A5007 /H/3 ~ 1 . 

There are two main reasons why our analysis yields differ- 
ent results from those published by Mason et al. (1996). First, 
KUG 1031+398 having a redshift of - 0.043, the [N n]A 6583 
line coincides with the atmospheric B band. When correcting 
for this absorption feature, the [Nil] true intensity is recov- 
ered (Fig. |l]b) and our red spectrum appears different from the 
published one; different line-ratios and widths are therefore not 
unexpected. 

Second, the line-profile analysis of Mason et al. differs 
from ours in that, while we force each Balmer component to 
be associated with forbidden lines having the same velocity 
and width, Mason et al. allow these parameters to have differ- 
ent values for the Balmer and forbidden line components. As 
a result, they found three H/3 components (a narrow, an inter- 
mediate and a broad one), as well as two [O III] components 
(a narrow and an intermediate one); they also detected three 
Ha components (again a narrow, an intermediate and a broad 
one), but only a single [Nil] component (narrow). The mea- 
sured width of the narrow H/3 component is 150 ± 20 kms -1 
FWHM, while the width of the narrow [O III] lines is 265 ± 10 
kms -1 ; this last value, significantly larger than the narrow H/3 
line width, suggests that the [O III] lines may have a complex 
profile. Moreover, the width of the [N II] lines is found to be 
significantly larger (400 ± 60 kms -1 ) than that of the narrow 
Ha component (190 ± 40 kms -1 ); this could be due to an in- 




Fig.2. Enlargement of part of the 
CCD frame showing the H/3 and 
[Om]AA4959, 5007 lines in the 
spectrum of KUG 1031+398 (aver- 
age of three 20 min exposures, af- 
ter cosmic rays removal); the ex- 
tended (narrow) component of H/3 
is clearly visible. 



accurate correction of the atmospheric B band, as we have seen 
above. 

Although our spectra have a lower resolution than those ob- 
tained by Mason et al. (3.4 A compared to 2 A FWHM), this 
does not affect the analysis; the narrow core components being 
identified and subtracted, all the discussion is centered on the 
broader components, well resolved even with our lower reso- 
lution. Similarly, the larger slit width used in our observations 
(2'.' 1 compared to l'/5 for Mason et al.) does not affect the study 
of these broader components, since only the contribution from 
the extended emitting region (the H II region, Fig. |2|), removed 
with the core, changes with the slit width. 

Boiler et al. (1996) and Wang et al. (1996) have suggested 
that the small width of the broad Balmer lines and the soft X- 
ray excess characteristic of NLS 1 galaxies could be the effect 
of a high accretion rate on an abnormally small mass black 
hole. Mason et al. (1996) have argued that, although the emis- 
sion line spectrum in KUG 103 1+398 is dominated by the ILR, 
a weak broad component is present with line-widths of the or- 
der of 2 500 km s -1 FWHM and that, therefore, at least in this 
object, such a model is not required. 

Our analysis of the spectra of KUG 1031+398 has shown 
that, in the BLR, the Balmer lines are well fitted by a 
Lorentzian profile with ~ 1 060 kms -1 FWHM; this value is 
much narrower than the value found by Mason et al. (~ 2500 
kms -1 ). This is due to the fact that we used a Lorentzian, 
rather than a Gaussian profile to fit the broad Balmer lines; 
the Lorentzian profile was required by the presence of broad 
wings, fitted with a Gaussian by Mason et al. (1996). 
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We have shown (Goncalves et al. in preparation) that in 
NLSls the broad component of the Balmer lines is generally 
better fitted by a Lorentzian than by a Gaussian; the Lorentzian 
Balmer lines (component 4), without any measurable associ- 
ated forbidden line, would qualify this object as a NLS1 with, 
in fact, very narrow lines. So, in this respect, KUG 1031+398 
is a normal NLS 1 and could be explained by the same small 
black hole mass model as the other objects of this class. 

Table 1. Emission line profile analysis of KUG 1031+398. 
I(H/3) and I(Ha) are in units of 10 -16 ergs _1 cm^ 2 . In 
columns 1 and 2 we give the mean of the relative velocities 
and widths measured on the blue and red spectra. The FWHMs 
are corrected for the instrumental broadening. 
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A6583 
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1.27 
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- 95 
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10.2 


0.84 
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81 


3 


- 395 


1115 


6.1 




18 




4 


- 105 


1060 






189 


938 



5. Conclusions 

We have obtained new spectra of KUG 1031+398 around H/3 
and Ha. We have found that the emission-line spectrum of 
this object can be modelled with four components: an extended 
H II region, two narrow emission regions of Seyfert 2-type and 
a relatively narrow "broad line" component, well fitted by a 
Lorentzian profile. 

We disagree with Mason et al. on the analysis of the emis- 
sion line profile of KUG 1031+398, in the sense that we find 
no evidence for the presence of an "intermediate" component in 
which the forbidden lines are almost, but not completely, sup- 
pressed by collisional de-excitation. Nevertheless, we find that 
this object is exceptional in having a "narrow" line region (de- 
fined as a region where A5007 /H/3 > 5) with almost the same 
width at half maximum as the "broad" line region (Balmer lines 
with no detectable associated forbidden lines); however, in the 
first case, the line-profile is Gaussian, while in the second case, 
it is Lorentzian. 
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